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We report entanglement swapping with time-bin entangled photon pairs, each constituted of a 
795 nm photon and a 1533 nm photon, that are created via spontaneous parametric down conversion 
in a non-linear crystal. After projecting the two 1533 nm photons onto a Bell state, entanglement 
between the two 795 nm photons is verified by means of quantum state tomography. As an important 
feature, the wavelength and bandwidth of the 795 nm photons is compatible with Tm:LiNb 03 -based 
quantum memories, making our experiment an important step towards the realization of a quantum 
repeater. 


INTRODUCTION 


Entanglement swapping entangles two photons that 
have no common past [l|. This fascinating phenomenon 
not only stimulated curiosity to understand quantum cor¬ 
relations but also plays an important role in vari¬ 

ous applications of quantum information science, includ¬ 
ing quantum computing [8. 9] and quantum repeaters 
[lQj . A quantum repeater-based communication chan¬ 
nel, for instance, exploits entanglement swapping to en¬ 
tangle interim nodes in a heralded fashion, and connect 
elementary entangled links (connecting interim nodes) to 
distribute ent ang lement in principle over arbitrarily long 
distances [ll|, [l2|. In turn, the resulting entanglement 
can be used to generate a secret key between distant users 

St 

Necessary ingredients for quantum repeaters, in addi¬ 
tion to entangled photon pairs and entanglement swap¬ 
ping, are optical quantum memories 0. Such memories 
allow the reversible mapping of (entangled) states be¬ 
tween light and atoms, and thereby remove the necessity 
for all elementary links to establish entanglement simul¬ 
taneously. While entanglement swapping has been re¬ 
ported before 00, the spectra of the resulting entan¬ 
gled photons were either orders of magnitude too large, 
or their wavelength were not suitable to allow subsequent 
interfacing with optical quantum memory. Here we re¬ 
move this impediment: using cavities, we spectrally en¬ 
gineer swapped photons in time-bin qubit states in such 
a way that their wavelengths (around 795 nm) and band- 
widths become compatible with our solid-state quantum 
memories 0- Furthermore, we accordingly increase the 
coherence times of the photons used to swap the entan¬ 
glement (both around 1533 nm wavelength), which will, 
in the future, allow such quantum interference measure¬ 
ments even with photons that have traveled through tens 
of kilometres of deployed standard telecommunication fi¬ 


bre [22j. To verify successful entanglement swapping, 
we measure (conditional) two-photon visibility curves. 
For more complete information, we additionally employ 
quantum state tomography to derive the concurrence and 
the fidelity of the swapped state with the nearest max¬ 
imally entangled state, and we compare the latter with 
the predictions of a recently developed model [23|. As 
with the use of quantum state tomography, this has not 
been done previously for time-bin qubits after entangle¬ 
ment swapping. 


EXPERIMENT 

A schematics of our experimental setup is depicted in 
Fig. □ A 1047 nm wavelength laser emits 6 ps long 
pulses at 80 MHz repetition rate. After second har¬ 
monic generation (SHG) in a periodically poled lithium 
niobate crystal (PPLN), the now 18 ps long pulses, cen¬ 
tered at 523.5 nm wavelength, travel through an unbal¬ 
anced Mach-Zehnder interferometer (MZI) whose path 
length difference corresponds to 1.4 ns travel time dif¬ 
ference, thereby splitting every pulse into two. Pairs of 
pulses emitted from the two outputs of the interferome¬ 
ter then pump two 10 mm-long PPLN crystals, in which 
spontaneous parametric down-conversion (SPDC) leads 
to time-bin entangled qubits |24) encoded into pairs of 
photons with wavelengths centred around 795 nm and 
1533 nm. Assuming, for the sake of explanation, for 
the moment that only individual photon pairs are cre¬ 
ated, this process yields states of the form \<&+)ab = 
(|e)^ 4 |e)_g + \£)a\£)b) emitted from one crystal, and 

|<3 >~)cd = -^{\e)c\e)D ~\£)c\2)d) emitted from the 
other crystal, where \e) and \t) represent early and late 
time-bin qubit states. 

Using a grating monochromator connected to a single 
photon detector, we measure the spectral width of the 
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FIG. 1. Schematics of of our setup. See text for details. 


795 nm photon at full width half maximum (FWHM) 
to be 1.5 nm. By energy conservation, and taking into 
account the bandwidth of the pump photons, we calcu¬ 
late the spectral width of the 1533 nm photons to be 
5.6 nm (FWHM). These values by far exceed the maxi¬ 
mum bandwidth of 10 GHz (corresponding to 21 pm at 
795 nm wavelength), over which quantum memories have 
so far been reported to operate [21|. Hence, to allow fu¬ 
ture interfacing with memories, we filter the 795 nm pho¬ 
tons using Farbry-Perot cavities (FP; one per source) to 
6 GHz, and also reduce the bandwidth of the 1533 nm 
photons to 12 GHz (again using Fabry-Perot cavities), 
which corresponds to 94 pm. (The filter cavities fea¬ 
ture only one spectral order within the spectral width 
of the created photons.) This additionally ensures that 
the photons’ coherence time, around 37 ps, exceeds that 
of the pump pulses, as required for entanglement swap¬ 
ping. The cavities are based on Corning ultra-low ex¬ 
pansion glass spacers and fused silica mirrors to reduce 
thermal noise, and their transmissions are 35% and 90% 
for 795 nm and 1533 nm photons, respectively. Using 
the approach described in ( 25 }, we find the probability 
of having a photon pair per source and pump laser pulse 
after filtering, i.e. per qubit, to be about 19.1 ± 1.8%. 

To swap entanglement to the two 795 nm qubits 
(propagating along spatial modes A and D), the two 
1533 nm qubits (in modes B and C) are subjected to 
a so-called Bell-state measurement (BSM) after travel¬ 
ling through short standard telecommunication fibres. 
This measurement is performed by sending the two pho¬ 
tons into the different input ports of a 50/50 beam 
splitter. Provided they exit through two different out¬ 
put ports and in different temporal modes (one pho¬ 
ton early, the other one late), their joint state is pro¬ 


jected onto the ) B c = 775 (|e)s|^)c - K)e|e)c) Bell 
state, leaving the two 795 nm photons in the entan- 
gled |^+}ad = ^(\e)A\^D + \e)A\e) D ) Bell state. We 
note that, in principle, it is also possible to make pro¬ 
jection measurements onto the |T + ) Bell state, thereby 
increasing the measurement’s efficiency from maximally 
25% to 50% To detect the 1533 nm photons, we 
use free-running, tungsten silicide (WSi)-based super¬ 
conducting nanowire single-photon detectors (SNSPD) 
cooled to around 0.9 Kelvin [I?}. However, due to fibre 
transmission loss inside our cryostat, we find a system 
efficiency of around 50%. Furthermore, we measure a 
detection-time jitter of 250 ps (FWHM), which is suffi¬ 
ciently small to allow resolving the temporal qubit modes 
(spaced by 1.4 ns), and the dark-count rate of 10 Hz en¬ 
sures little noise-pollution of detection signals. 

For the Bell-state measurement to work, the two 
1532 nm photons must be indistinguishable at the beam 
splitter, i.e, their spatial, temporal, spectral, and polar¬ 
ization modes must be identical. This is verified using 
so-called Hong-Ou-Mandel (HOM) interference (28|: if 
two indistinguishable photons (not qubits) impinge on 
a symmetric beam splitter from different input ports, 
then they bunch and leave together by the same output 
port due to destructive interference between the prob¬ 
ability amplitudes associated with both input photons 
being transmitted or both reflected. Conversely, if the 
two photons are distinguishable, no such interference oc¬ 
curs, and they leave the beamsplitter with 50% proba¬ 
bility through different outputs, resulting in coincident 
detections. The HOM visibility, defined as Vhom = 
(A "max A min ) /A max |29| , where A max crrM A min denote 
coincidence count rates measured with completely dis¬ 
tinguishable and (maximally) indistinguishable photons, 
respectively, is a common way to characterize the degree 
of indistinguishability. We find Vhom = 27.5 db 2.5%. 
This value is consistent with the fact that our 1533 nm 
inputs are not single photons, but are mixtures of pho¬ 
ton Fock states with thermal distribution, for which the 
HOM visibility is upper bounded by 1/3 [29j. 

Repeating the measurement conditioned on the de¬ 
tection of two 795 nm photons (described below) and 
with small pump power, which projects the 1533 nm 
inputs onto approximate single photons, the HOM vis¬ 
ibility increases up to 87.5 ± 5.5%. We attribute the 
gap to the theoretical value of 100% to insufficient spec¬ 
tral filtering of the photons, and, to a much smaller 
extent, to remaining contributions from the simultane¬ 
ous emission of more than two photons. Indeed, using 

Vhom = l//l + (^) 2 ® EH , where AT denotes the 
duration of the pump pulses, and r the coherence time of 
the 1533 nm photons after filtering, we find that the max¬ 
imum visibility achievable with our setup is 89%, which 
corresponds to our result within experimental error. 

To assess the quality of the entangled state after the 
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entanglement swapping, the two 795 nm photons are 
sent through Mach-Zehnder interferometers that intro¬ 
duce the same travel-time difference of 1.4 ns as the in¬ 
terferometer that acts on the pump beam and thus al¬ 
low projecting photon states onto various time-bin qubit 
states [24[. To ensure phase stability during the mea¬ 
surements, i.e. constant projectors, the interferometers 
are passively temperature stabilized. Additionally, their 
phases are actively locked using a frequency-stabilized 
laser at 1550 nm wavelength that is also sent through 
the interferometers, and a home-made feedback loop. Fi¬ 
nally, the 795 nm photons are detected using four stan¬ 
dard silicon avalanche photodiode-based single photon 
detectors with efficiencies around 50%, detection jitter 
of 500 ps, and dark counts around a few hundred Hz. 
All detection signals are recorded using a time-to-digital 
converter that is started by a successful BSM. 

We also measure the total heralding efficiency (32| . 
which, for the 795 nm (1533 nm) photon, we define as the 
ratio between the rate of photon pair coincidence detec¬ 
tions and the rate of the detection of 1533 nm (795 nm) 
photons. Hence, for the low mean photon number in 
our experiments, the heralding efficiency is equal to the 
probability of detecting a photon that has been created 
by downconversion in the source. This takes into ac¬ 
count all transmission loss (such as loss in filter cavities, 
prisms and fibre coupling loss), non-unity detection effi¬ 
ciency, and a fundamental restriction caused by the band¬ 
width mismatch between the pump and downconverted 
photons. A simple calculation of the bandwidth restric¬ 
tion (see appendix) shows that the latter alone limits the 
heralding efficiency for the 795 nm (1533 nm) photons 
to 17.4% (34.8%). Considering furthermore the known 
values of 50% (70%) for detection efficiency, 40% (80%) 
for transmission through the cavity, and 85% (85%) for 
transmission through the prism, we anticipate a herald¬ 
ing efficiency of 2.96% (16.6%) for the 795 nm (1533 nm) 
photons. Experimentally we obtain heralding efficien¬ 
cies of 1.96% for 795 nm and 5.8% for 1533 nm pho¬ 
tons. The difference between the measured and expected 
values is predominantly due to imperfect fibre coupling, 
which we assess to be 66% and 35% for the 795 nm and 
1533 nm photons, respectively. Hence, we find that one 
of the main limitations to achieving higher heralding ef¬ 
ficiencies is the bandwidth mismatch between the pump 
photon and the downconverted photons. This can be im¬ 
proved by using a spectrally narrower pump, as further 
described in the Appendix. 

MEASUREMENTS AND RESULTS 

To investigate the presence of entanglement between 
the two 795 nm photons after the BSM with the 1533 nm 
photons, first, we measure the probability for coincidence 
detection of two 795 nm photons in detectors A and D, 
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FIG. 2. color online) Coincidence detection probabilities 
(conditioned on a successful Bell-state measurement) of the 
795 nm photons as a function of the phase difference in the 
measurement interferometers. The uncertainty bars indicate 
one standard deviation and are calculated from measured 
rates assuming Poissonian detection statistics. The visibility 
of the sinusoidal fit is 56.2 ± 5.7%. The blue curve features a 
visibility of 1/3, indicating the classical bound for a separable 
Werner state. 


conditioned on projecting the two 1533 nm photons onto 
the |T _ ) Bell state, as a function of the phase difference 
a — /3 of the two analyzing interferometers. This cor¬ 
responds to projecting each 795 nm photon onto equal 
superpositions of \e) and \t). (For these as well as the 
following measurements, the mean photon pair number 
was set to 19.1%.) The results, depicted in Fig. [2j show a 
sinusoidal curve with a fitted visibility of (56.2 =b 5.7)%. 
The phase difference axis is recalibrated based on the fit 
so that a — /3 = 0 corresponds to a maximum of the co¬ 
incidence probability. The sinus period is fixed to 2i r, 
however, refitting the data with the period as a free pa¬ 
rameter yields an almost identical curve with a period 
of (0.98 =b 0.03)27r. This confirms the excellent stabil¬ 
ity and calibration of the relative interferometer phases 
during our measurements, which required a total of 36 
hours (the average four-fold coincidence rate was about 
10/hour). Furthermore, the visibility clearly exceeds the 
maximum visibility of 33% that can be obtained using a 
separable Werner state [33|. This confirms the presence 
of entanglement, provided the often made assumption of 
having a Werner state is satisfied. 

To remove this assumption, we reconstruct the den¬ 
sity matrix of the 795 nm photons’ joint quantum state 
by means of maximum likelihood quantum state tomog¬ 
raphy (QST) Q, which will allows assessing additional 
measures that quantify entanglement and derive more 
information about experimental imperfections. Towards 
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this end we perform a total of 36 joint projection mea¬ 
surements corresponding to all combinations of projec¬ 
tions onto eigenstates of <rx, cry and cr^. Each combina¬ 
tion requires a coincidence measurement (conditioned on 
a successful Bell-state measurement) with specific (local) 
interferometer phase settings (for projections onto eigen¬ 
states of ax and ay), or measurements of photon ar¬ 
rival times (for projections onto eigenstates of oz\ We 
emphasize the importance of interferometer phase sta¬ 
bility, which we verified above, to ensure proper (and 
stable) measurements. The reconstructed density ma¬ 
trix is depicted in Fig. [3J It allows calculating the con¬ 
currence (C), an entanglement measure that is zero for 
a separable state and larger than zero for an entangled 
state [35j. We find C = 0.36 =b 0.07, which exceeds the 
threshold by five standard deviations (for this and the 
further experimental values the uncertainty is calculated 
by means of Monte Carlo simulation and assuming the 
coincidence counts follow Poissonian statistics), thereby 
confirming the conclusion of having an entangled state 
derived from the visibility in Fig. O We also compute 
our reconstructed state’s fidelity with the expected state 
|T + ) to be F = (68 db 3)%. Optimization over all max¬ 
imally entangled states also shows that |T + ) yields the 
largest fidelity with our reconstructed state (within the 
step size of the optimization algorithm). 



FIG. 3. color online) Density matrix characterizing the joint 
state of the two 795 nm photons after entanglement swapping. 
The fidelity with respect to the |T + ) Bell state - which we find 
to be the nearest maximally entangled state - is (68 ±3)%. 

Furthermore, we can now assess how well the recon¬ 
structed density matrix is described by a Werner state 
by calculating the fidelity of the swapped state with the 
nearest Werner state a = + \{1 — v)l. Here, 


|V’) is any maximally entangled two-qubit state, 1 is the 
identity matrix and v is a variable that parametrizes the 
state. Evaluating F = (tryj\fap\fa) 2 over all \x/j) and 
f, we find that the maximal value F = 97.5+4 4 % is 
achieved for \ijj) = |T + ) and v = 0.592. Hence, using the 
fidelity measure indicates that the reconstructed state is 
close to a Werner state. The remaining difference is due 
to the admixture of non-white noise in the reconstructed 
state, which we primarily attribute to the residual dis- 
tinguishability in the Bell-state measurement (this im¬ 
perfection only affects measurements of a x and a y - not 
a z [HI)- Hence, it is reasonable to use the fitted visi¬ 
bility in Fig. [2j which is based on measurements in the 
a x and a y bases, to infer the presence of entanglement 
in the state after entanglement swapping. However, the 
observation that our state is not an exact Werner state 
also leads to a note of caution, in that a visibility mea¬ 
surement comprising projections onto a z - which in our 
experiment is less affected by experimental imperfections 
- may provide an inflated measure of the entanglement 
in the swapped state. 

To independently assess how experimental imperfec¬ 
tions limit the amount of entanglement in the final bi¬ 
photon state, we employ the method described in [23] 
with measured experimental parameters for heralding ef¬ 
ficiency (1.96% at 795 nm and 5.8% at 1533 nm), HOM 
visibility (89%), mean photon pair number (19%) and fi¬ 
delities (95%) for the individual sources. This leads to a 
density matrix having a concurrence of 0.43, i.e. a state 
that is slightly closer to a maximally entangled state than 
the measured one. We find the fidelity of the density 
matrix predicted by this model and the one we measured 
to be 98.1+4.0%, which suggests that we understand the 
noise sources affecting the swapped state. Because this 
model accounts for non-white noise arising from the lim¬ 
ited HOM visibility in the BSM, the predicted density 
matrix features a slightly larger fidelity with our mea¬ 
sured state as that resulting from the comparison with 
the closest Werner state. 

Based on our analysis, we identify multi photon pair 
emissions, due to the probabilistic nature of SPDC 
sources, as our main source of errors that limits the value 
of the concurrence. To improve our result, we therefore 
have to decrease the probability of generating a photon 
pair per pulse below our current value of 19%. For in¬ 
stance, reducing the mean photon pair emission probabil¬ 
ity to 10% would, according to our model [23], result in 
an increase of the concurrence to 0.53 (keeping all other 
parameters unchanged). Note that alternative way to re¬ 
duce the impact of multi-pair emissions are to exploit a 
quantum Zeno blockade to suppresses multi photons [37[ 
or a quantum non-demolition measurement that reveals 
the number of simultaneously emitted photon pairs [HI, 
thereby allowing in theory to completely ignore detec¬ 
tions stemming from multi pair emissions. According to 
the model this would yield a concurrence of 0.82. 
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OUTLOOK AND CONCLUSION 

Before concluding, let us briefly discuss the possibil¬ 
ity to extend the current implementation into an ele¬ 
mentary quantum repeater link, in which case the two 
1533 nm photons have to travel tens of km before being 
submitted to the Bell-state measurement, and the two 
795 nm photons have to be stored in optical quantum 
memories. We note that these criteria are easily met as, 
first, their long coherence time makes the 1533 nm pho¬ 
tons robust against travel-time fluctuations during long¬ 
distance transmission, second, loss in optical fibres at this 
wavelength is minimal, and third, as quantum memories 
for 795 nm photons with 6 GHz bandwidth are avail¬ 
able 0. However, even assuming memory efficiencies 
exceeding 50%, which remains to be demonstrated for 
memories of such large bandwidth, the coincidence count 
rates are currently too small to demonstrate an elemen¬ 
tary quantum repeater link - let alone building a useful 
one. Solutions, on the one hand, are relaxed focussing of 
the 523 nm laser pulses into the SPDC crystals, which 
has been shown to significantly improve the coupling ef¬ 
ficiency ( 39 ]. Secondly, the lossy cavities for the 795 nm 
photons can be removed as the memories themselves will 
work as filters. The latter will increase the coincidence 
rates by one order of magnitude. Furthermore, given the 
large bandwidth of the SPDC photons (of which we cur¬ 
rently use only a few GHz), and the large bandwidth and 
spectral multi-mode storage capacity of, e.g., atomic fre¬ 
quency comb-based Tm:LiNb 03 quantum memories, it 
is possibility to work with many frequency channels in 
parallel using a quantum repeater architecture that we 
described in [12j. 

To summarize, we have experimentally demonstrated 
the creation of entanglement between two 795 nm pho¬ 
tons, whose properties are suitable for further storage in 
broadband quantum memories such as our Tm:LiNb 03 
waveguides, by means of a Bell-state measurement 
with two 1533 nm photons, each of which was initially 
entangled with one of the two 795 nm photons. Our 
demonstration constitutes an important step towards 
the generation of a quantum repeater: the heralded 
entanglement of two quantum memory-compatible 
photons. 
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APPENDIX: HERALDING EFFICIENCY 
VIS-A-VIS PHOTON BANDWIDTH 


The bandwidth mismatch between the pump light and 
the filters used in our experiment for the downconverted 
photons fundamentally limits the maximum attainable 
heralding efficiency for each of the down-converted pho¬ 
tons (commonly referred as signal and idler). In this 
section we develop a simple model - inspired by the pic¬ 
torial representation in [40| - that captures the main con¬ 
sequences of this effect. 

As described in the main text, we define the heralding 
efficiency for the signal (tjh s ) as the ratio between the 
coincidence detection rate (C%) and the single detection 
rate for the idler (Si): 


Vh s 



(i) 


Similarly, the idler heralding efficiency is found as: 


Vh, = 



( 2 ) 


The pump laser has a spectral bandwidth Av p given 
by the pump pulse duration. Since the down-conversion 
process must conserve energy, the frequencies of any pair 
of down-converted photons are correlated to within the 
pump bandwidth. As a result, the Joint Spectral Am¬ 
plitude (JSA) of the down-converted photons (see Fig. 
E} can be illustrated by a diagonal band (green region) 
with a cross sectional width corresponding to the spec¬ 
tral width of the pump (Av p ). Note that we assume 
that, over the relevant bandwidth, the phase-matching 
condition is less restrictive than that resulting from en¬ 
ergy conservation, and thus we do not indicate it on the 
figure. However, phase-matching determines the total ex¬ 
tent of the diagonal band (outside the view of the figure) 
and thus sets the overall bandwidth of down-converted 
photons. 

Next we explore the effect of restricting the bandwidth 
of either the signal or the idler photon by means of spec¬ 
tral filters with bandwidths A v s and Az/$, respectively. 
In Fig. m and for the case in which we filter the idler, 
this corresponds to carving out a horizontal band in the 
JSA. The idler photon single detection rate is propor¬ 
tional to the overlap of this horizontal band with the di¬ 
agonal band representing energy conservation (red area 
in Fig. 0^), i.e. Ai = \/2Au p Aui. Similarly, for the case 
of filtering the signal, the single count rate is propor¬ 
tional to the area of the blue region in Fig. 0)3, which 
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is A s = \[2Av v Av s . These expressions are approximate 
because they assume that the count rates of the filtered 
photons are directly proportional to the specified areas. 
This, however, is only correct if the spectral profile of 
the filters and the JSA given by the pump bandwidth 
are box-shaped (flat-top). Moreover, the area is only cal¬ 
culated as above if the filter bandwidths are smaller than 
the pump-bandwidth i.e. Ais s w < A v v . 

Finally, we extend this picture to the case in which 
we filter both signal and idler and measure the coinci¬ 
dence count rate. As illustrated by Fig. the coinci¬ 
dence count rate is proportional to the area (in pink), 
i.e. the intersection of the horizontal and vertical bands 
given by the filters A s i = Av s Ai Clearly, this area is 
smaller than both A s and A^ thus limiting the heralding 
efficiency. More precisely the heralding efficiencies can 
be expressed as: 

_ Csi _ A s i _ Av s 

VHs ~ 

Csi Asi Av{ /r%\ 

m = sr = X = V2ZT, (3) 

In our experiments Av s = 6 GHz (795 nm photon), 

Av s = 12 GHz (1533 nm photon) and Av p = 24.4 GHz 
(523 nm photon). Using these values in Eq. [3j we obtain 
t]h s = 17.4% and r\u i = 34.8%. Though these values are 
calculated using a simple model, they provide a quali¬ 
tative explanation of the maximum attainable heralding 
efficiency for each of the signal and idler photons in our 
system and allows us to assess the loss due to other fac¬ 
tors such as optical elements, optical fibre coupling, and 
detector efficiency. 
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